function [xpct_mean] = pctile_mean(x,pct)
 x=sort(x);
 pct=ceil(pct*size(x,1));
 xpct=x(pct);
 xpct_mean = [mean(x);xpct];
end